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Abstract 

Making use of symmetry properties and known exact results for 
the Hubbard Model on bipartite lattice, we show that (1) there is 
no phase separation for repulsive coupling at low dopings, (2) phase 
separation and superconductivity co-exist in the ground states for at- 
tractive coupling in a range of filling fractions. 
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Recent experiments indicate the presence of phase separation in cuprate 
superconductors, namely, there exist separate hole-rich and hole-poor phases |l|] 
[^]. There have been theoretical studies, all employing approximations or 
numerical estimations, on this phenomenon based on the t-Jmodel0-|p, as 
well as directly on the Hubbard model Most of these studies indicate 
that these models can explain the phenomenon of phase separation. The 
exceptions are those of Putikka et a/. 0, and Moreo et al. [|^. While the 
quantum Monte Carlo calculations of small 2D lattice systems indicate no 
phase separation for the 2D Hubbard model 0, the investigation using high- 
temperature expansions for the 2D J model suggests the existence of phase 
separation only for a limited range of J/t values]^. 

In this note, we look again at this phenomenon on the basis of the 
Hubbard model with repulsive on-site interaction. We make use of some 
exact results for the Hubbard model on bipartite lattice with attractive 
coupling |^-[|10|, which, after applying the well-known unitary particle-hole 



transformation p^, imply a certain general property for the repulsive cou- 
pling case. It in turn implies that the repulsive Hubbard model on bipartite 
lattice can not accommodate phase separation at low dopings. We also look 
at the case of attractive Hubbard model, for which we show that, at tem- 
perature T = 0, phase separation is possible, and actually co-exists with 
sup er conduct ivity. 

The Hamiltonian of the Hubbard model on bipartite lattice is 

ij,a i i 

with the usual notations. The bipartite lattice sites A are divided into subsets 



A and B. The number of sites of (B) is denoted by 1^41 and the total 

number of sites by |A| . In (1), tij vanishes if both i and j belong to the 
same subset. An exact result obtained by Lieb]^ states that, for attractive 
coupling {U < 0) and even number of electrons, the ground state of the 
system is unique and has total spin S = 0. Thus, at temperature T = 0, 
the system is non-ferromagnetic. This result was extended to T > by 
Kubo and Kishi|Ty], who find that there is an upper bound on a certain two- 



point function. It implies absence of long-range magnetic order in attractive 
Hubbard model at finite temperature. 

We introduce an external magnetic field in the 2-direction. The Hamilto- 
nian (1) becomes 



Define 



s' - If, (3) 

where sf is the z-component of the spin of the particle on site i. In general, 
we have 

hm {s^)=f{h,T), (4) 

|A|— >oo 

where (...) is the average over a grand canonical ensemble. Invoking Lieb- 
Kubo-Kishi's result, we have 



/(o,r) = o, (5) 

and, furthermore, since there is no long-range magnetic order, that f{h,T) 
is analytic at /i = 0. As a consequence, we can draw the conclusion that 



f{h,T) is a well-behaved, i. e. continuous and single-valued, function in the 
neighborhood of h = for T > 0. For small h, 

lim {s^) = x{T)h, (6) 

where x{T) is the spin susceptibility. 



We now make the unitary particle-hole transformation [|n| 

Ci| Qt, Cii eicl^, (7) 

where = 1 (—1) on the lattice set A (B). The Hamiltonian (2), up to a 
constant, is transformed into 

(ij),'^ « (8) 

U' = -U, 

^-f^' = h, (9) 
/.' = f -/X, 



with 



and (6) into 



hm D = l + x(T)(f/'-2/i'), (10) 

A — »oo 



where D is the density for the system described by H, 

d=?^^4t^- (11) 

|A| 

Recall that U' = -U > 0. 

From (10), we see that for a system described by a repulsive Hubbard 
model {U' > 0) on bipartite lattice, the density D in the thermodynamic 



limit is a well-behaved, i. e. continuous and single-valued, function of the 
chemical potential near namely at low dopings. This behavior is in 
contrast with what is expected of a first-order phase transition. 

When a system with chemical potential fio undergoes a first-order phase 
transition at temperature T, and separates into two phases with different 
densities, we have 

Di = lim^^^„+o limiM^oo D, 

[12) 

D2 = lim^^^„_olim|A|-.oo^, 
with Di 7^ D2. The density D, as a function of the chemical potential, is 
discontinuous at /io in the thermodynamic limit when there is first-order 
phase transition. Thus the continuity of lim|A|^oo-D in (10) implies that a 
repulsive Hubbard model on bipartite lattice can not accommodate co-existence 
of two phases at low dopings. This result is valid for general h' , including 
h' = 0, which is the case of vanishing external magnetic field. 
The following are two remarks: 

1. The conclusion we made with regard to the phase separation in a re- 
pulsive Hubbard model on bipartite lattice applies to low dopings. This is 
because the exact result of Lieb-Kubo-Kishi is for the case of vanishing mag- 
netic field [h = 0). However, if the function f{h,T) is also well behaved, as 
may be reasonably expected for a paramagnetic system, at large values of h, 
then our conclusion about phase separation can be extended to high dopings. 

2. Lieb has proved in the case of repulsive Hubbard model on bipartite 
lattice with half-filling and h = that the (degenerate) ground states have 
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the same spin S = ^{\A\ - \B\) (if \A\ > \B\), with 

S. = -1{\A\ - \B\), -\B\) + 1, - \B\) (13) 

This result can be used to make certain statement concerning phase sepa- 
ration at T = for attractive Hubbard model on bipartite lattice. Making 
use again of the unitary particle-hole transformation, it is straightforward to 

obtain, for the attractive case, that at T = 0, 

\A\ - \B\ AA\ - 151 1 \A\ - \B\ 

This implies that in an attractive Hubbard model on bipartite lattice different 
phases can co-exist at T = 0, with the densities having values given in (14). 

It has been shown by Shen and Qiu|l2| that for attractive Hubbard model 
on bipartite lattice with |v4| — |S| = 0(|A|), the ground states with densities in 
the range (1 — ^^j^^) < D < (1 + ■'^^j^^) have off-diagonal long-range order 
(ODLRO) in the two-particle reduced density matrix, and, consequently, 
posses the property of superconductivity [|13|. That is, superconductivity and 
phase separation co-exist in this case. 

In conclusion, we have investigated the question of phase separation for 
the Hubbard model on bipartite lattice. Our result that there is no phase 
separation at low dopings for T > 0, being based on the exact results of 
Lieb0, and Kubo-Kishi||lOl, is an exact statement of the repulsive Hubbard 



model on bipartite lattice of any dimension. It is consistent with the quan- 
tum Monte Carlo calculations of small 2D lattice systems for the Hubbard 
Model0]. It is also consistent with the high-temperature expansion result of 
Putikka et a/. for the 2D t-J model, to the extent that the Hubbard model 
and the t-J model are expected to have similar results at small J/t values. 
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